Glassy dynamics and hysteresis in a linear system of orientable hard rods 
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We study the dynamics of a one-dimensional fluid of orientable hard rectangles with a non- 
coarse-grained microscopic mechanism of facilitation. The length occupied by a rectangle depends 
on its orientation, which is coupled to an external field. The equilibrium properties of our model 
are essentially those of the Tonks gas, but at high densities, the orientational degrees of freedom 
become effectively frozen due to jamming. This is a simple analytically tractable model of glassy 
phase. Under a cyclic variation of the pressure, hysteresis is observed. Following a pressure quench, 
the orientational persistence exhibits a two-stage decay characteristic of glassy systems. 

PACS numbers: 



I. INTRODUCTION 

Almost all liquids, if cooled sufficiently fast, form a 
glassy structure; much effort has been directed towards 
understanding this phenomenon. Monte Carlo simula- 
tions have shown that systems with purely hard-core 
interactions can describe qualitatively much of the ob- 
served phenomenology of the glass transition, e.g., the 
very fast rise of relaxation times, and the absence of an 
associated latent heat. Theoretical analysis of such mod- 
els is hampered by an incomplete understanding of the 
equilibrium (fluid-solid) phase transition. The best stud- 
ied hard-core model is a system of hard spheres, which 
has experimental realizations in colloidal, granular and 
other systems [H-Q- I R many cases size dispersion, or 
other built-in complexity, is introduced in order to avoid 
crystallization and thereby observe a glass transition 0- 
0]. For monodisperse systems, there is no transition in 
one dimension [8[, while two and three dimensional sys- 
tems are highly prone to crystallize. In higher dimen- 
sions nucleation rates are low and the glass state 
is more easily attained. Indeed, in the limit of very high 
spatial dimensionality, there seems to be an ideal glass 
transition [(| [l3| , although how far this extends down to 
lower dimensions is still debated uM. 
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Hard core potentials need not be spherically symmet- 
ric. On a lattice, for example, where the symmetry is dis- 
crete, the behavior strongly depends on the dimension of 
the system, the lattice structure, and the exclusion range 
(see [la, 13 and refs. therein). 

Here we study the dynamic properties of a one- 
dimensional system of classical hard rectangles, with 
only two orientations allowed, horizontal and vertical. 
These will be called "rods" in what follows. Classical 
linear fluids have been extensively studied over the last 
decades d, [r7rl2e| ; the case of elongated rods with ori- 
entational freedom has also been analyzed [2(J [29|-[34j]. 
While the equilibrium properties of this class of models 
are well understood, their dynamic properties, in partic- 
ular those related to the glass transition, have received 
somewhat limited attention [35M39| . The simple, one- 
dimensional model considered here, reproduces (at least 
in part) the glassy phenomenology with an explicit (non- 
coarse-grained) microscopic mechanism of facilitation. In 
analogy with closely related models in which the con- 
straints are kinetic instead of geometric [4fj| - |42j . due to 
a spatial and temporal coarse-graining, here there is no 
thermodynamic transition and the only nontrivial behav- 
ior is kinetic. In the present context this represents an 
advantage, in that the relaxation is not complicated by 
critical slowing down or metastability. 

At short time scales, the glassy phase can be modelled 
as a metastable phase, in restricted thermal equilibrium. 
An exact calculation of the properties of such metastable 
states, including the equation of state, has been made 
recently for a toy model [3^, [43|. In these papers, the 
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crgodicity is explicitly broken, and one assumes that in 
the glassy phase some transition rates are exactly zero. 
The simple model considered here provides an extension 
of the treatment in Ref. [39| to include the description of 
slow evolution of the macroscopic glassy state at longer 
time scales. In our model, ergodicity is not explicitly 
broken, and the dynamics is capable of bringing the sys- 
tem to equilibrium, but at high pressures, the relaxation 
of orientations becomes so slow that these variables are 
effectively frozen over any reasonable time scale. Un- 
der steady increase of pressure, they never fully relax. 
Macroscopic properties in the frozen regime are found to 
depend on history, in particular, on the rate at which the 
pressure is increased. 

The structure of the paper is as follows. The next sec- 
tion introduces the model while Sec. IIIII reviews its equi- 
librium properties. We describe the dynamic behavior in 
section HVl Sec. [V] contains some concluding remarks. 



II. MODEL 



We consider a system of N rigid rods of length a > 1 
and unit width on a line of length L. Each rod is de- 
scribed by the variables (a;,, Si), where Xi € [0,L) is the 
position of the center of mass and Si denotes its orienta- 
tion (0 for horizontal, 1 for vertical). Rods in state Si = 
occupy a length a and those in state Si = 1 occupy a unit 
length (see Fig. [I]). Since the order of the particles cannot 
change, we take x\ < X2 < ■ ■ ■ < x^. For convenience, 
we define Xq = and xn + i = L; there are no orienta- 
tional degrees of freedom associated with these variables. 
We define inter-rod distances, = Xi + \ — x^, so that 
Si=o Vi = L. When all rods are constrained to have 
the same orientation or, cquivalcntly, a = 1, the origi- 
nal Tonks gas 0] is recovered. Without this restriction, 
the model is analogous [44j to a binary mixture of rods 
(with conservation of the total particle number, but not 
the number of each species separately) . The particles are 
subject to an external field h' , coupled to the orientations 
{Si}, so that the potential energy takes the form 



n = J2<t>bc(yi) + h'J2 s i 



(1) 



where the first term denotes hard core interactions and 
the sums extend over the N rods; for h' > the horizon- 
tal orientation is favored. 

The hard-core interaction between the ith rod and its 
right neighbor is 



rods i and i + 1, given by 



1 



a+-{Si+S i+1 )(l-a). 



(2) 



An important quantity is the mean free volume per par- 
ticle, Vf = v — a + (a — l)m, where v = p^ 1 = L/N and 
m is the fraction of rods in the vertical orientation, or 
"magnetization" . 
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FIG. 1: Example of a rod configuration. 

We define a local, continuous-time, stochastic dynam- 
ics for the positions (or equivalently, the separations yi), 
and orientations. The time evolution is Markovian, and 
satisfies the detailed balance condition. Each rod exe- 
cutes unbiased diffusive motion. The wall at xn+i = L 
also undergoes diffusive motion, but is subject to a bias, 
with ratio of jump rate that increase L by SL to the back- 
ward transition is cxp(— pSL). The rods can also change 
their orientation, rotating about its center of mass. While 
the orientation varies continuously during a transition, 
we suppose the latter occurs so rapidly that orientations 
other than vertical or horizontal may be neglected. In 
addition, the dynamics includes an important constraint 
on orientational transitions: particle i can only change its 
orientation if its neighbors to the left and right are suf- 
ficiently far away that its rotation (through 90°) is not 
blocked geometrically (see Fig. [2]). In order for particle i 
to change its orientation, it is necessary that 



Vi > 



S i+ x(l-a) 



(3) 



where r = s/l + a' 2 is the rod diagonal. The analogous 
relation for yi-\ must be satisfied as well. Note that 
these conditions depend on the states of the neighboring 
rods (i ± I), but not on S% itself. Without the nonover- 
lapping constraint, r is replaced by a in the above equa- 
tion. Any violation of the excluded volume condition is 
rejected. These transitions are accepted in accordance 
with the Metropolis criterion, that is, with probability 
min[l,exp(-/3AH)]. One Monte Carlo step (MCS) con- 
sists in an attempt to update all degrees of freedom (i.e., 
all positions and orientations, and the volume). 

Our Monte Carlo simulations were performed on sys- 
tems of N — 1000 rods; and results represent averages 
over 100 (or more) independent realizations. 



0hc(j/i) 



0, y% > a l 



^oo, yi < Oj, 

where a; is the minimum distance between the centers of 



III. EQUILIBRIUM PROPERTIES 

The equilibrium properties of the model are indeed 
very simple. The canonical configurational partition 
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FIG. 2: Geometrical constraint associated with an orienta- 
tional transition. Analogous restrictions apply also for the 
other possible orientations. 



Eqs. (JT)) and (|5J) imply the relation 

v = - + m eq + <t(1 - m eq ), (9) 

which implies that the free volume per particle is Vf = 
l/p, as in the Tonks gas. It is worth noting than in 
equilibrium, the variables {y^} and {Si} are all mutually 
independent. The behavior of m as a function of the 
pressure for several values of a and h is illustrated in 
Fig. El 



function is given by 

Z N (h,L) = Ve-^* 5 ' ! dyx---dy N (4) 
{St} Jr 

where h = f3h' and the subscript R denotes the restric- 
tions Y^iLiVi = L an d Vi — ^ a ii where the minimum 
distances a* are defined in Eq. @. The usual factor of 
1/N\ is absent due to the fixed order of the particles on 
the line. Introducing the variables z% = yi — Oj, (with 
zq = x\ and Zn = L — xn), we have, 



Z N (h,L) = Y^e-' l ^ s * [dz --- [ 
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dz 



N 



xS 



N 



N 



J2zi-L + Ncr-(<T-l)J2Si 



i=0 



1=1 



We study the system in the constant-pressure ensemble; 
the partition function is 

Y N (h,p) - J dLe-P L Z N (h,L), 

where p denotes the pressure divided by fc^T. A simple 
calculation yields 



Y N (h,p) = e- pNa p- N (l + n) 



(5) 



with k = exp[(cr — l)p — h]. The Gibbs free energy per 
particle is given by g = — (A r /3) _1 InF, where, in the 
thermodynamic limit, 



lim — In Ym 



-op + ln(l + k) — hip. 



In this limit the volume per particle is, 
1 dlnY 



v(p, h) 



-!)«,! 

) 

P 



N dp 1 + k 

while the fraction of vertical rods is, 

. , , 1 <91nF k 

m eq (p,h) 



N dh 



1 + K 



(6) 



(7) 



(8) 




FIG. 3: Equilibrium fraction of vertical rods as a function 
of pressure p, for h — and 4, (empty and filled symbols, 
respectively). Points are simulation results, while the lines 
correspond to Eq. ([8]). When p is large enough, all rods are 
aligned vertically. 



IV. TIME-DEPENDENT PROPERTIES 

In this section we study the kinetics of the magneti- 
zation m and molecular volume v. We assume that the 
diffusive relaxation is much faster than the orientational 
relaxation. Then the displacement degrees of freedom 
may be assumed to be in thermal equilibrium. At high 
pressures, most update attempts for orientation change 
fail, as they are blocked. Suppose first that the orienta- 
tions are fixed. In this case, the different sectors in the 
pico-canonical ensemble are specified by orientation of 
each rod. Different sectors with the same number of verti- 
cal rods are macroscopically equivalent. Within a sector, 
the displacement degrees of freedom are assumed to be in 
equilibrium. Then the translational dynamics will bring 
the system to a constrained equilibrium distribution, in 
which the distances zi are independent, exponentially dis- 
tributed random variables with mean l/p, and the mean 
length of the system is (L) = N{l + (l-m)a} + (N+l)/p], 
where m is the fraction of vertically oriented particles, 
not necessarily equal to the equilibrium value m eq . 

Now, allowing the orientations to fluctuate, an equa- 
tion of motion for m(t) can be derived if we assume that 
the translation dynamics is rapid, so that between any 
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pair of successive oricntational transitions, the interpar- 
ticle distances attain the constrained equilibrium distri- 
bution mentioned above. Under this hypothesis, the par- 
ticle orientations are mutually independent, as in equilib- 
rium (we assume as well that the orientations are initially 
uncorrelated) . 

Consider, for example, a transition from Si = to 1. 
This would be allowed only if the two gaps on the two 
sides of the rod are large enough. In the restricted equi- 
librium ensemble, separations are independent, exponen- 
tially distributed random variables, P(y) = pcxjp(-py), 
and the probability that a hole larger than w = (r — a)/2 
appears at one of its sides is P(vf > w) = cxp[— p(r — 
f)/2]. Given that such gaps must exist at both sides and 
that this flip is against the field, we see that the effec- 
tive transition rate is proportional to exp[— p(r — a) — h], 
and that it is independent of the state of the neighboring 
rods. Thus we may write the transition rate for Si = 
to 1 as 



7+ = jcxp[-p(r -a) - h], 



(10) 



where 7 is an arbitrary attempt rate, independent of p 
and h. On the other hand, if the transition is from 1 to 
0,w = (r- l)/2 and 



7- 



7 exp[— p{r — 1)] 



(11) 



Note that only 7+ depends on h and that these rates 
satisfy detailed balance. For small values of the pressure, 
free space is abundant and the slowest process is a flip 
against the field, so that the larger time scale involved 
is given by 7+. On the other hand, when the pressure 
is large enough, the production of large enough holes is 
the dominant slow process and the relevant characteristic 
time now scales as jZ ■ 

The evolution of m(t) is governed by 

dm 

— =7+(l -m) -7_m= -rm + 7 +! (12) 

where r = 7+ + 7_. If the rates are time independent, 
then letting <f> = m — m eq =771 — 7+ /I\ we have 



A particularly interesting example is that of a system 
initially in equilibrium at pressure po, and subject to a 
pressure that increases linearly with time, p{t) = po + Xt 
for t > 0, where A is the annealing rate. In this case, 



G(t) 

so that 

lira G(t) 



,—h—po(r—a) 



1 



,-\(r-<7)t 



+e -po(r-l) 



r — a 
!_ e -A(r-l)t 



,—h—po(r—a) 



1 



j-poO'-l) 

r - 1 



(16) 



(17) 



which is finite for A > 0. Inserting this result in Eq. (|14p . 
we see that the initial magnetization mo is not "forgot- 
ten" even when t — > 00. It is easily verified that memory 
of the initial magnetization persists for a pressure in- 
crease of the form p(t) = po + Xt a , for any positive values 
of A and a. 




FIG. 

5, 



Fraction m of vertical rods versus pressure p for h = 
a = 2, 7 = 1 and several values of the rate of pressure 
increase, A. From top to bottom, A = (equilibrium, bold 
line), 10 -3 , 10 -2 , and 10 -1 . Points: simulation; solid lines: 
theory, Eq. (fl4)) . 



dt 



■r<t>, 



(13) 



showing an exponential approach to equilibrium. 

The system is driven out of equilibrium if the pressure 
(or the external field, h) is time-dependent. Suppose that 
7 + and 7_ depend on time through the pressure. Then 
we have 

m(t) = e- g{t) (m + J ds 7+(s) e e < s ^ , (14) 



where 



g(t) = f dt'T(t'). 

Jo 



(15) 



Examples of m(t) (for pressure increasing linearly with 
time) are shown in Fig. 2] for Po = 1, h = 5, a = 2, and 
7 = 1. For A = 10~ 4 , the difference between m{p) and 
the equilibrium result is small. For larger rates of pres- 
sure increase, on the other hand, there are marked differ- 
ences between the final value of m and the equilibrium re- 
sult. Although m(t) is well described by Eq. (jl~4")) . as can 
be seen by the excellent agreement with the simulation, 
the same does not occur for the molecular volume, v(p), 
and the free volume, at larger rates of pressure increase 
A. The reason is that in this case the rate of pressure in- 
crease is large enough so that we can no longer treat the 
interparticlc distances as in instantaneous thermal equi- 
librium, on the time scale of the oricntational relaxation. 
This effect can be incorporated in an approximate man- 
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FIG. 5: Volume per molecule v versus pressure p for the same 
parameters as in Fig. U during the pressure annealing. Points 
represent the simulation data while the solid lines are the 
improved theoretic predictions based on Eq. (| 18f) (with 7 ~ 
0.3). The lines for A = 10~ 4 and A = (equilibrium, bold 
line) are indistinguishable at this scale. Inset: free volume 
vf- Analogous deviations from equilibrium are again seen for 
large values of A. 



ner in the theory if we assume that the free volume per 
molecule follows a relaxational dynamics, that is, 

£~ r -(—;) < 18 » 

To evaluate the transition rates 7+ and 7_ we require 
the probability density p(z). The simplest hypothesis is 
that p(z) is exponential, as in equilibrium, but with the 
mean free volume in place of its equilibrium value, 
1/p, so that p(z) = Vp 1 exp(—z/vF)- Then the evolution 
of the magnetization is given by Eq. (|12[) , with transition 
rates as in Eqs. (|TU|) and (jTTj) . but with p replaced by 
1/vp. With an appropriate choice of the relaxation rate 
r„, this simple theory yields reasonable agreement with 
simulation results at larger quench rates, as is shown in 
Fig. [5] Some deviations between the theoretical predic- 
tion and simulations are evident at the highest quench 
rate, for smaller pressures; this is not surprising given 
the simplifications introduced. 

In addition to the "freezing out" of the orientational 
degrees of freedom under a steady pressure quench, the 
system exhibits interesting hysteresis effects under an os- 
cillatory pressure. Hysteresis loops obtained through nu- 
merical simulation are illustrated in Fig. [BJ for various 
values of A = \dp/dt\ in triangle- wave cycles of pressure 
variation. Similar results are obtained via numerical in- 
tegration of Eq. (|12[) . assuming rapid equilibration of the 
free volume. Notice that for larger rates, the system de- 
scribes a sequence of irreversible loops before entering a 
reversible one. For even larger rates than those shown in 
the figure, we find a greater number of irreversible loops, 
analogous to those obtained in compaction experiments 
of rods under vibration 1451 . 



FIG. 6: Fraction of vertical rods under cyclic variation of 
the pressure, showing hysteresis loops. For the highest rate, 
A = 0.1, the loop only closes after several cycles. Parameters: 
h — 5, a — 2; the pressure varies between 1 and 20 at rates A 
as indicated. 



It is worth noting that an external field h is not re- 
quired to observe freezing. Even with h = 0, the rapid 
reduction in the transition rate 7+, Eq. (|10p . with in- 
creasing pressure ensures that the orientational degrees 
of freedom cannot equilibrate. Inhibition of orientational 
relaxation is greatest for \/l + a 2 —a as large as possible, 
i.e., for a tending to unity. (Of course this tends to re- 
duce the excess of v over its equilibrium value.) Thus, for 
h = 0, a = 1.1, and other parameters as above, one finds 
moo = hmt-j.00 m(t) = 0.620 if the pressure increases at 
a rate of A = 1, and ?7ioo = 0.5409 for A = 10. (As A 
is increased, TOqq approaches the initial magnetization, 
equal to 1/2 for h = 0.) 

If instead of a smooth annealing, the system is sud- 
denly quenched from low to high pressure, a two-step 
decay typical of glassy behavior is observed in the ori- 
entational persistence (the fraction of rods that did not 
flip since t = 0). Fig. [7] shows the results for several val- 
ues of the final pressure after the system is quenched 
from an equilibrium state at po = 1. Whatever the 
value of the final pressure, the initial decay of the per- 
sistence function is exponential and ruled by the against 
the field flip of rods that are horizontal at t = 0, so that 
P(t) ~ exp(— 7+i) in this regime. As the final pressure 
is increased, the flip from 1 to becomes slower and 
a plateau develops, after the initial fast decay, whose 
width increases with the final pressure. The height of 
this plateau is on the order of the initial magnetization, 
too- (Its precise value is slightly smaller than mo since 
some of the up rods will have already flipped during the 
initial fast regime). When the final pressure is large, it 
takes a long time for a rod to flip, since enough space 
must be freed, which requires a cooperative rearrange- 
ment of the neighboring rods. This happens with a rate 
proportional to 7_. Thus, the data in Fig. [7] are well fit 
by a sum of exponentials, corresponding to the fast and 
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FIG. 7: Persistence P(t) as a function of time as the system 
pressure is quenched from 1 to p at t = 0, for cr = 2 and 
h — 5. When p is small, the decay is exponential while for 
larger values, after an initial exponential decay, the persis- 
tence develops a plateau whose height is close to the initial 
magnetization mo and whose width increases with p. The 
solid lines are from Eq. (|19l) . 

slow processes in the model: 

P(t) ~ e" 7 +' +m e- fc7 - t (19) 

where, for the parameters of Fig. fc ~ 0.4 and mo ~ 
0.01. The coefficient of the second term, 0.01, is the 
height of the plateau and roughly corresponds to the ini- 
tial magnetization. Thus, starting with an even smaller 
initial pressure, and thus a larger magnetization, the 
plateau can be tuned to higher values. The width of the 
plateau increases with the final value of the pressure and 
diverges as p — > oo. This diverging relaxation time, 7Z 1 , 
is related to the increasing length of the cooperative re- 
gion [33j]. Notice that, at variance with other models for 
the glass transition, the slow relaxation is not associated 
with stretched exponentials. 

V. CONCLUSIONS 

One-dimensional systems with short-range interac- 
tions, such as the model studied here, do not exhibit 



a phase transition at finite temperature and pressure in 
equilibrium. A jamming transition invol ving certain de- 
grees of freedom may nevertheless occur [33|, with a di- 
verging length scale, as the control parameter (inverse 
pressure or temperature) goes to zero. Here we study a 
simple, geometric model on the line, for which analyti- 
cal results for both the statics and the dynamics may be 
obtained and compared with numerical simulations, with 
excellent agreement. This system of hard rods is subject 
to geometric constraints that prohibit a rod changing its 
orientation if the distance from its neighbors is too small 
(i.e., in the absence of sufficient free volume). In our 
model, the time spent in transit between the two orien- 
tations is assumed to be small, and the orientation degree 
was taken to be a discrete variable. It is straightforward 
to extend the discussion to the case where we allow con- 
tinuous orientations, although no qualitative differences 
are expected. 

In the presence of an external field that disfavors the 
vertical position, our model may be seen as an in-layer 
description of a higher-dimensional system. In a dense 
system, the vertical position will be disfavored due to 
excluded- volume interactions with the rods in the neigh- 
boring layers. Thus the external field may be interpreted 
as an effective interaction to take into account the re- 
maining dimensions. 

Despite its simplicity, the model presents several prop- 
erties characteristic of the glass transition, such as an- 
nealing rate dependence and two-step relaxation. The 
dynamical behavior can be understood in terms of the 
two microscopic reorientation processes involved. 
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